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^ . Abstract 

$-H ' 

O , 1 The hydrodynaniic equation derived by N-particle statistical mechanics is investigated. The reason 

.^*H ' is an attempt to provide additional information concerning the closure problem of turbulence theory. 

, In this paper,the third-order term of the friction force of the equation is calculated in terms of third- 

00 ■ and fourth-order equilibrium time correlation functions. 

T-H ■ 

'Th '' 1 Introduction 

^' . . 

This is the first of a series of papers which deal with the hydrodynaniic equation derived by N-particle sta- 
tistical mechanics; especially with the calculation of the third order term of the friction force (divergence 

5][)' of the stress tensor). 

jyT ' The motive for the investigation is the closure problem of classical turbulence theory. It is well known 

O ■ that the problem arises when one declares, in case of turbulent fluid motion, the hydrodynaniic velocity 

' jy^ ' u to be a random variable, and attempts to obtain an equation for its expectation by just building 

^~>. the expectation over the Navier-Stokes equation. The non-linear term then gives rise to an additional 

1-5 ' variable; when one formulates an additional equation for this, another variable is obtained, an so on. 

Apparently, if the basic equation of motion is non-linear, this method does not work. Instead, from the 
point of probability theory, it would be necessary to find the probability density /(u, t) of the turbulent 
process, and to calculate the expectation as: 

OO 

m- (u)= / u/(u,i)du (1.1) 



o 



Ph' 



o 

• I / ist a functionally dependent on u; thus, (11.11) is a functional integral. As early as 1952, Hopf has 

'^ ■ formulated this theoretical approach [Ij; a short account can be found in [5], ch. 4. - It became apparent 

that it is extremely difficult to provide an equation for (u) by following this approach. On the other 
hand, the success of direct numerical simulation (see, e. g. [3j) leads to the conclusion that it must be 
possible to calculate the mean velocitiy on the basis of the Navier-Stokes equation. 

In this situation, it is perhaps reasonable to attempt the derivation of the mean velocity equation by 
^^ I N-particle statistical mechanics, in order to provide further information on the nature of the problem. It 

5h ■ is obvious that for investigation of turbulent processes it is not necessary to deal with the molecular na- 

. . .' ture of the fluid. It is the probabilistic structure of statistical mechanics which constitutes the difference 

to Navier-Stokes theory. Statistical mechanics is the basis for describing a large set of physical processes 
including fluid dynamics; it would be a surprise if turbulent motion would be an exception. The compli- 
cation is that in a statistical mechanics analysis it is necessary to distinguish between macroscopic and 
microscopic parts of the motion, and to handle the latter in a suitable way. - The Navier-Stokes equation 
can be derived from statistical mechanics via the Boltzmann equation by Chapman-Enskog method; see 
e. g. [4]. But in the course of the derivation of the Boltzmann equation the molecular chaos assumption is 
introduced, which is correct for many applications including laminar fluid motion, but seems problematic 
for a turbulence investigation. 

An essential step is defining the mean velocity in the frame of statistical mechanics. The fiuid velocity 
is defined: 

u=i(p) (1.2) 

P 



p is the mass density, p the microscopic momentum density. In this paper, we consider incompressible 
constant density and temperature processes. Then u is, up to a constant factor, equal to the expectation 
of p. Under rather general conditions, the expectation is equal to the arithmetic mean of a time scries 
of the microscopic quantity. Since for laminar flow fluctuations are microscopic, u then is the "point" 
velocity of the flow as we usually understand it. On the other hand, in turbulent flow there appear 
macroscopic fluctuations; the arithmetic mean by definition averages over these also; thus, in this case u 
is already the mean velocity. We have the peculiar situation that the statistical mechanics u is equal the 
classical hydrodynamic u or (u), depending on whether the flow is laminar or turbulent. Therefore, the 
hydrodynamic equation in the form derived by N-particle statistical mechanics will, by way of trial, be 
considered the mean velocity equation. 

The derivation by means of Zwanzig-Mori projection operator technique of statistical mechanics 
(POT) is utilized, in the presentation of Grabert [S]. For the convenience of the reader, some introductory 
material from an earlier paper of the author [5] is repeated. The hydrodynamic equation obtained by this 
technique shows a formula for the friction force containing a local equilibrium time correlation, which 
is a non-linear functional of the velocity; presently, there exists no theory from which it could be calcu- 
lated. Therefore, it is necessary to develop it into a functional power series in the velocity; this presently 
restricts the applicability of the equation to low Reynolds number flow, just outside the Navier-Stokes 
regime. The coefficients of the series now contain total equlibrium correlation functions, which have been 
calculated by employing a separate theory. 

In this paper, calculation of the third order term of the power series is described. Applying correlation 
theory is left to a consecutive paper. Performing a test with empirical data is still under way and will be 
reported later. 



2 Hydrodynamic equation 

This section is a somewhat altered version of section 2 of [6j. There, the description was restricted to 
stationary processes. In this paper, I provide the usual time-dependent formulas. The fluid is considered 
to be a system of N particles of mass m with positions y^ and velocities Vj which are combined to the 
phase space matrix z. Vector components are described by Latin indices, e. g. Yj — {yja\ ■ The system 
is enclosed in a box of Volume V . Later on, the thermodynamic limit is performed. - A function g{z) 
is called a phase space function, or microscopic variable. Especially, we need the space densities of the 
conserved quantities mass, energy and momentum n, e, p which are collected to a 5-element linear matrix 
a. They are functions of an additional space variable x: 

N 

"■ = '^'^3^i^-yj) (2-1) 

For the particle functions Uj we have rij — m , pj — m'Vj , while the energy function contains the 
interparticle potential, m ist the particle mass. Quantities (like a) which are lists of 5 elements are 
denoted by normal letters while three-component vectors and tensors are bold. - The quantities a obey 
the conservation relations: 

a = -V-s (2.2) 

The fiuxes s have the same general structure as the a (j2.ip : especially, we have Si = p . The time 
evolution of any phase space function A is described by the Liouville equation: 

A ^CA (2.3) 

C (often defined as iC) is the Liouville operator. From (|2.3p . the formal solution for A{t) given the 
initial value A is: 

A{t) = e^* A (2.4) 

In the statistical model, z and N are considered random variables; that is, the probability density 
f{z,N) is of grand canonical type. The ensemble mean value (expectation) of a phase space function A 
is defined in the 'Heisenberg' picture: 

{A)it) = E / dzA{z,N,t)f{z,N) (2.5) 



In this formula, f{z, N) is the initial probability distribution, and A{z, N, t) is the value of A at time 
t if the initial positions and velocities of the particles are described by z . The operation (integration + 
Summation) is sometimes indicated by the symbol 'tr': 

tr{f7} = 51 / dzn{z,N) (2.6) 

Certain probability densities (also called distributions here) are frequently used in the analysis. One 
of them is the (total) equilibrium distribution which corresponds to macroscopic rest: 

/o = 4iN) exp($o + Kf^N - H{z))) (2.7a) 

m) = ^^i^r"" (2.7b) 

Here, h is Planck's constant, /? = l/{kBT) , ks being Boltzmann's constant and T the temperature, 
fi is the chemical potential which is a function of mass density p = (n) and temperature, and H{z) 
is Hamilton's function which describes the total energy of the fluid. For the normalization constant, 
we have $o = ~(3Py , P being the equilibrium pressure. Expectations with respect to the equilibrium 
distribution are denoted by ()o. In case of a simple fluid, the 'relevant probability distribution' of Grabert's 
formalism (see [5], sec. 2.2) is the local equilibrium distribution: 

fL{t) = 4>{N) exp($(i) - a{z) o h{t)), (2.8a) 

6={/5(^«'--),/3,-/3u}, (2.8b) 

2 m 

$(f) = -log(tr{V'exp(-ao6(t))}). (2.8c) 

Here the symbol o is introduced for the operation: Multiplication, plus Summation over the 5 elements 
of the linear matrices a , h, plus Integration over geometrical space. The elements of h are called the 
conjugate parameters; they are functions of the quantities /3 , p, and u which we will sometimes call the 
thermodynamic parameters, and which will be considered to be slowly varying functions of space and 
time. The h are deflned such that the expectations of the a are identical to their expectations in local 
equilibrium: 

{a) = (a)i (2.9) 

The POT is a means for separating macroscopic and microscopic parts of the random variables. It 
starts by deflning the set of phase space functions which are relevant for the description of the process. 
For simple fluids, this set is identified with the densities of conserved variables, a . A projection operator 
is deflned which projects out of any microscopic variable A the part which is proportional to the relevant 
variables. It reads: 

VA = {A)l + {A5a)L o {5a 5a)l^ o 5a (2.10) 

Here, 5a — a — (a) l', ()2 denotes the inverse of the expectation matrix in the formula. For general 
non-stationary flow, the local equilibrium distribution (j2.8ap . and therefore V, are time-dependent. Then, 
instead of [S] (2.11), for the decomposition of e-'-* the formula [S] (2.4.1) is obtained: 



e£t ^ e^* V{t)+ / dt' e^*' r{t'){C - P{t')){l - P{t'))g{t' , t) + {I - V{t))g{0, t) (2.11) 

Jo 

g{t',t) = exp_( / dt"C{l-V{t"))) (2.12) 

Jt' 

Q{t' ,t) is a time-ordered exponential Operator which describes the time dependence of the dissipative 
part of the equation of motion. The analysis in [5] consists in applying (|2.4[) to a, using (|2.11[) . By aver- 
aging over the initial probability density, and after some manipulations, Grabert's generalized transport 
equation 5], (2.5.17) ist obtained. It is postulated in POT that the initial probability density is of the 
form of the 'relevant probability density', which for simple fluids is defined to be the local equilibrium 
density (j2.8ap . Grabert states that this should not be considered a general restriction of the method but 
a means to form the general particle system into the type specially considered (the simple fiuid here); 
see [5], sec. 2.2 . It is shown that in this case the last term in (|2.1ip vanishes after averaging. - Finally 



the part of the formula pertaining to the momentum density is taken, and (|1.2p is used. The result, [5] 
(8.1.12), (8.1.13), reads with the donotations employed here: 

p(^+u- Vu) = -VP + V-R (2.13) 



R(x,t)= / dt' /"dx'S(x,x',i,t')V'u(x',t') 



(2.14) 



S(x, x', t, t') = pmt', i)(l - 7'(i)s(x)](l - P(t'))s(x'))L,t' (2.15) 

The microscopic Definition of the pressure P{x,t) [5j (8.4.7) will not be repeated here. {)L,t denotes 
an expectation with respect to the local equilibrium distribution at time t. Equation (|2.13|) has the formal 
structure of the hydrodynamical equation; but at the present state of the analysis, it is still an exact 
equation. It is the approximations performed later on in the formula for the stress tensor R which will 
transform it into an irreversible equation. - In Addition to (|2.13p . the equation for the particle number 
density component of (a) yields the continuity equation: 

V-u = (2.16) 

The kernel function S of the stress Tensor i? is a time correlation function in local equilibrium which 
is a nonlinear functional of u. As far as the author knows, presently it is possible to calculate correlation 
functions for total equilibrium only. Therefore, it has been necessary to expand S into a functional power 
series in u. As has been explained in section 3 of [6], the expansion can be done with the set b{t) of 
conjugate parameters of the local equilibrium formula (|2.8b[) . The expansion is performed at the point 
b = bo, which corresponds to u = 0: 

6o = {-/?-, /3,0} (2.17) 

m 



.1 .2 



b^bo = {(3-u\0,-l3u} (2.18) 



The power expansion of S reads: 



S = S|,„ + Jk*(6-6o) + ^^ko**{(^'-^'o),(&-6o)} + ••• 
= S(°) + S(i) + S(2) + ... (2.19) 

The designation of terms in the second row is for later reference. The * indicates multiplication, 
summation over five elements and integration over space and time. For the present purpose, the expansion 
is cut after S*-^-'. When these terms are inserted into (J2.14L one obtains corresponding parts R*-^-*, R*-^', 
R('^^ of R; the index again describing the order in u. The elements of bo are the conjugate parameters 
of the total equilibrium ensemble. Thus, when S and its derivatives are taken at 6 = 6o, the quantities 
in the integrand resemble total equilibrium space-time correlation functions. In [B], these quantities have 
been calculated by Kawasaki technique [7]; though this theory still contains certain intuitive elements. 
In the present state of the project, for correlation functions the multilinear mode-coupling theory (MCT) 
of Schofield and co-workers has been used [S]. For the 3-point correlations in S'^', the results of the two 
methodes coincede; for 4-point quantities, there are differences. - In order to calculate S^^', it is sufficient 
to restrict b — bo to the last element of (|2.18p . since the first element is of second order in u; but for S'^', 
this element of S^^^ has to be added. 

It is an important test for the calculation method that (I2.13p . correctly approximated, should yield the 
Navier-Stokes equation. The linear part R^"'^^ of the stress tensor has been calculated by several authors 
including Grabert [S], sec. 4; the Stokes form of the stress tensor is obtained, with a microscopic definition 
of the friction matrix. There remains another detail of the argumentation: The left-hand side of equation 
(j2.13p is of second order in the velocity. Therefore, in order for the equation accurately derived, it must 
be shown that the second-order part R^^-* of the stress tensor vanishes. The present author published 
several papers in ArXiv (S] dealing with this topic, and finally could show that R*^^) = 0, so that the 
Navier-Stokes equation derives correctly. 



3 Third-order term of the stress tensor 

Derivation of the third-order term provides considerable effort; a computer algebra system (Mathematica) 
has been used for most of the calculations. In order to describe the third-order term, we switch from the 
stress tensor R to the friction force D = V • R. By introducing the corresponding term of (|2.19|) into 
(j2.14p . we obtain a formula for the third-order term D3 of D. Again, only the u-part of (|2.18p has to be 
considered. For the detailed investigation, we need the indexed form of the formula: 

(D3)a(x,i) = ^V, f'dt' fdx' r dt" fdx" r dt'" fd^"' ^^^-bcd{^^<t,t') ^ 

^ ^ ^ ' 2 'Jo J Jo J Jo J SKM",t")Sbf{jc'",t'")^'" 

X V>6(x',i')ue(x",i")u/(x"',i"') (3.1) 

Latin letters denote indices which run over 3 elements; by contrast, greek indices run from 1 to 5. - In 
order to continue, the second-order functional derivative of S is to be calculated. By |S] (second paper) , 
the first-order derivative consists of two parts. The parametric part is given by (A7), already specialized 
to u = 0; here we need the general form: 



SSa: 



bed 



,(x",i") 



= -Sit" ~ t')(3 {[Gif ,t)Q{t)saci^)]Q{t')Sa,ix" ,t')Qit')sMi^'))L,t' (3.2) 

Q{t) = l-P{t) (3.3) 



Sa^{-K,t) = ae{-x.) - {a^{-x.))L,t (3.4) 

The functional part is [9j (second paper) (A14): 

SSabcd 



(x",i") 



^(3e{t"-t')e{t-t"){[g{t\ncr{nsa,{^",nQ{ng{t'',t)Q{t)sac{^)]Q{t')sM{^'))LM 

(3.5) 
For later reference, the first-order derivatives are formulated with respect to all five conjugated param- 
eters (the greek index e is used). - These formulas, as well as the starting formula (I2.15|) . contain operator 
chains which are, as usual, written like products. The question arises whether I can, when differentiating 
a chain, use the product rule; in other words, if A and B are operators, whether it is correct to write: 

SAB SA^ .SB ,^^, 

Since I could not find a general proof in the literature, I checked p.6p for any two consecutive operators 
appearing in p.2p . (j3.5p . and found it correct in all cases. The calculations are partly long and tedious and 
will not be shown here. - From p.2p . (13. 5p . the second-order functional derivatives are to be calculated. 
The formulas depend on 6 by 6 and 8 'factors', respectively; including, in each formula, the dependence 
of the local equilibrium probability density. The corresponding terms are given in Appendix A, (jA.2[) to 
(jA.lsp . In the course of the calculation, for CJ{t',t) an identity formula is needed which can be checked 
by discretizing it: 

Q{t')g{t',t) = Q{t')g{t\t)Q{t) (3.7) 

Three of the terms are found to vanish; thus, the second-order derivative of S, and therefore the 
third-order term of R consist of 11 terms. Finally to be added are parts of the first-order derivative (|3.2p . 
p.Sp which stem from e = 1. - For application in (13. IL all terms have to be taken at 6 = bo- This has 
effects listed below: 

- expectations in local equilibrium change into those in total equilibrium, denoted by ()o; 

- projection operators V{t), Q{t) change to time-independent operators V, Q; 

- 5ae{x,t) changes to ae{x); 

- g{t',t") changes to the exponential operator g{t) with t — t" — t': 

git) = e^S* (3.8) 



This operator has to be substituted by the similar operator g{t): 



g{t) = e2^* (3.9) 

The reason is that for g{t) there is an identity which substitutes it further to the non-projected 
exponential operator. - The substitution succeeds for aU 11 terms with the aid of one of the following 
identities: 

Qgit) - 9{t)Q (3.10a) 

g(t)C = Cg{t) (3.10b) 

5(i)Q/: = -VC + Cgit) (3.10c) 

Next, we change from the microscopic variables a, s to the corresponding ortho-normalized quantities 
/i, r; this extracts from each term the factor {%Y ■ Together with the factor /3 which appears explicitly 

in each factor, the pre-factor in p.ip changes to ^p^P- The gradient operators appearing in (13.11) are 
transferred to the kernel function; by partial integration, this alters the sign of the expression. In 
the formula, we then have the expressions Vc'f'aci,'^), V^rM(x') which we again denote rQ(x), r(,(x') 
respectively; the distinction to the original r's is by the number of indices appended. Finally, we change 
to Fourier space. The expression now reads: 

-1 pt POO I^OO 

{D:,)^{t)^^-p^P dt' dt" dt"'K{ro,rlh*^,h;,t,t\t",nu,{t')u2{t")us{n (3.11) 

^ Jo JQ Jo 

Here we have introduced number indices. A number index is a combined index which contains an 
index and a wave number variable. Moreover, double indices contain a summation over the indexed 
elements and an integration over the corresponding wave number space, together with the factor ^^^,,3 . 
The kernel function K consists of the 11 Terms presented in Appendix A, (IA.16P to (|A.26I) . For shortness, 
I used the abbreviation f =Qr. 

4 Preparation of the kernel function 

Starting with this section all expectations are with respect to total equilibrium. We can therefore omit 
the at the expectation symbol. - The kernel function has to be brought into a form so that it contains 
time correlation functions which can be evaluated by MCT. These are of the general form {[f{t)A]B*)o 
mit f{t) = e^*. As a first step, the projection operators P and Q are to be substituted. As can be seen 
in the formulas in the appendix, the terms contain to 4 of these operators (the operator contained in f 
remains unaltered). With the reformulations done in the last section, the operators now read: 

rA={A) + {Ahl)hi (4.1) 

QA = A-{A)-{Ahl)h^ (4.2) 

Generally, the indices contained in the number indices in these formulas must be greek (running over 
the full list of 5 elements); it then would follow that I have to insert 6h4 instead of ^,4. Starting with 
the consecutive paper, we will restrict the calculation to shear mode flow (as opposed to sound and heat 
modes); it will then be seen that the indices in (j4.ip . (J4.2I) are latin. - Each substitution involves a new 
doubled number index. If several substitutions are to be done in the same term, the number indices have 
to be denoted differently. Therefore, the substitutions were performed in 4 consecutive steps. After each 
step, bracket products were expanded, and certain auxiliary formulas applied: 

{hif2) = {{Qhi)r2) - 

(CA) = 

{Cg{t)h) = (4.3) 

(r) = -{Ch) = 

Chi = — fi — (rift,4)/i4 

The third and fourth formula are special cases of the second. - The result for the 11 terms, now 
denoted Xi, • • • , Kn, is presented in the appendix, (jA.27|) . I continue to have the kernel function divided 
in this manner, for this makes checkups of the calculation somewhat easier. The 6- and ^-factors are 
temporarily omitted; they will be added again when time integration is considered. 



It is seen that the time correlations in these formulas belong to one of two groups described by: 

{[g{t)fo]B*) (4.4a) 

{[Cgit)ho]B*) (4.4b) 

i? is a certain microscopic function. The formulas are still built with the operator g{t) and are to be 
transformed into those with f{t). This will be done in a number of steps. We start with (|4.4ap . First, 
we have: 

{[git)ro]hl) = (4.5) 

p.lOap and Qhi — have been used. Next, we introduce a symbol for the memory function of the 
process: 

roi(t) # {[g{t)fo]rt) (4.6) 

The memory function is, as usual, "localized" in time, i. e. approximated by: 

roi(t)=7oi^(i) (4.7) 

7oi = {2TT)'S{ko - ki)7oi (4.8) 

7q]^, when inserted into (J3.11L for incompressible fluid reduces to (written with indices a, b): 

7oi == k^^ab (4.9) 

Sab is the Kronecker symbol. - Next, we need an operator identity described in [TU], Appendix B2; in 
the time-dependent form, and applied to the operators defined in this paper, it reads: 

g{t) = f{t) - [ dt'f{t')VCg{t - t') (4.10) 

Jo 

We obtain the correlation function identity, for any microscopic function B: 

{[g{t)fo]B*) = {[f{t)ro]B*) - [ dt'Toiit - t'){[f{t')hi]B*) (4.11) 

Jo 

After localization of the memory function, we get: 

{[g{t)fo]B*) = ([/(t)ro]S*) - i^Af{t)hA]B*) (4.12) 

The first term on the rhs of (j4.1ip is still to be transformed. The fifth formula (j4.3p together with 



(|2.4p and 0)04 = {rohX) yields: 

^0 = -dtho - W04/14 (4-13) 

It is necessary to write the constant matrices oj, k with a tilde since they show a form similar to (|4.8p . 
- We obtain: 

{[g{t)fo]B*) = -dt{[f{t)ho]B*) - k„4{[f{t)hi]B*) (4.14) 

Ko4=Wo4 + 7o4 (4-15) 

We denote with C„ die time correlation functions which are to be calculated by MCT: 

{CnU-n-l{t) = {[f{t)ho]hl ■ ■ ■ K_^) (4.16) 

Then, we have the following applications of (I4.14p : 

{[g{t)fo]hlh1,) - -dt{C^U2(t) - ^08(^^3)812(0 (4.17) 

{[g{t)fo]h\hlhl) = -9^(^4)0123(0 - ^08(^4)8123(0 (4.18) 

Somewhat more complicated is the transformation if B contains f as a factor. Again, (I4.13P is applied, 
together with uj^^ = — W04: 



{WYoVlK) = -dt{[f{t)h„]fih*) - ko4{[f{t)hi]fih;) 

= dt{[f{t)homhi)h;) - Q,,dt{[f{t)ho]h;h;) + ~Koi{{[f{t)hi]{dthi)h;) - uu{[f{t)hi]h;h;)) (4.i9) 

The first and the third term need further transformation. Formula (j3.1ip shows certain symmetries 
concerning time and wave number integration. Thus, the following formula is correct if it is only used 
within ([XTT]) : 

{dthl)h* = ^{{dthl)h* + hldth*) = ^dtiKh*) (4.20) 

We obtain, for the correlation function (j4.19p : 

= -^att(C3)012(i)-^175t(C3)072W-K04(i9t(C3)412(i)+^17(C3)472W) (4.21) 

We now switch to time correlations belonging to group (|4.4bp . We need an operator identity similar 
to (j4.10p which again can be found in [10], Appendix 2B. Applied to this paper, it reads: 

g{t) - fit) - f dt'g{t')V£f{t - t') (4.22) 

Jo 

We obtain a formula similar to (|4.1ip : 

{[£g{t)ho]B*) = {[Cf{t)ho]B*) - / dt'{[Cf{t - t')ho]hl){[£g{t')h4B*) (4.23) 

Jo 

For shortness, we introduce some new denotations: 

^o{t) = {[Cg{t)ho]B*) (4.24) 

$o(<) = {[Cf{t)h„]B*) = dt{[.f{t)ho]B*) (4.25) 

Then, applying (|4l^ . K23^ reads: 

*o(i) = *o(t) - / dt'dt{C2)o4{t - t')*4(t') (4.26) 

Jo 

This is an integral equation for ^ which can be formally solved by temporarily switching to Laplace 
space. The result is: 

*o(i) - Mt) + «04 / dt'$4(i') (4.27) 

Jo 

Taking (j4.25p into account, we can integrate (|4.23p : 

{[Cg{t)ho]B*) = dt{[f{t)ho]B*) + ko4{{[f{t)h4]B*) - (h^B*)) (4.28) 

This is the basic formula for time correlations of form (j4.4bl) . Taking B = hi, we have: 

{[Cg{t)ho]hl) = dtiC2)oi{t) + Ko7{{C2)7i{t) ~ Sn) (4.29) 

The lowest-order non-projected time correlation function C2 obeys a simple differential equation ([5], 
third paper, (1.11)): 

dt{C2)oiit) = ~ko7{C2)7iit) (4.30) 

This is introduced into (|4.29p : 

{[Cgit)ho]hl) = -^01 (4.31) 

The lowest-order time correlation of type (I4.4bp is a constant. - For the next-higher order, we find 
from (|T^ : 



{[Cg{t)ho]h\hl) = dt{C^U2{t) + So7((C3)712(t) - J712) (4.32) 

joi2 = (C3)oi2(i-0) (4.33) 

Again the situation is somewhat more complicated if B contains the projected flux f. For the lowest 
order, we obtain, by applying (14.28^ : 

{[Cg{t)h^]flh1,) = dt{[f{t)h^]flh1,) + hoMf{t)h4nh*^) - (h^rlh*)) (4.34) 

We apply (gH]): 

{[Cg{t)ho]fih;) = -dt{[f{t)ho]{dthi)h;)+uirdt{[f{t)ho]h;h;) + 

+koi{~{[f{t)hi]{dthl)h*) + Cj^r{[f{t)K]h*,hl) - {h^flh;}) (4.35) 

We introduce a notation for the static correlation appearing in the formula: 

Si24 = {fih2hl) = -{hh2hl)* = -{hirlK) (4.36) 

Finally, we can apply a symmetry consideration as in (14. 19^ : by using (|4.20p we obtain: 

{[Cg{t)ho]flh*^) = ^du{C3)oi2it) + CdudtiC3)o72{t) + 

+ko4{^dtiC3)4i2{t)+Cdn{C3)mit) + h24) (4.37) 

For the next-higher order correlation we have, by (14. 28^ : 

{[£g{t)ho]flh;hl) = dt{[f{t)ho]flh*hl) + ko4{{[f{t)h4nh;hl) - {h4rlhlhl)) (4.38) 

Introduction of (|4.13p yields: 

{[Cg{t)ho]rlhlhl) = -dt{[f(t)h^]{dthl)hlhl)+uudt{[f{t)h^]h*7hlhl)) 

+ka4{-{[f{t)h4]{dthi)h;hi) + uu{[f{t)h4]h;h;hi) - {h4fih;h;)) (4.39) 

The symmetry consideration which is correct if it is applied within the integral formula (|3.1ip gives: 

{dthi)h;hi = l{{dthi)h;hi + h\{dthi)hi + hih*^dthi) = \dt{h\hihi) (4.4o) 

Again a denotation for the static correlation is introduced: 

S4i234 = {fih2h3hl) = -{fih2h3hlY = -{h4flhlhl) (4.41) 

The final form for (|438| is: 

{[Cg{t)h^]flhlhl) = \du{C4U2z{t)+Co^7dt{C4)oi2:i{t) + 

+ ^04(|at(C4)4123(i) + Cjn{C4)472z{t) + S4i234) (4.42) 

The formulas obtained in this section are introduced into the terms (|A.27p . The Dirac and Heaviside 
functions omitted are added again. The kernel function K in p. lip now consists of 5 main parts: 

K = 5{t"' -t')5{t" -t')Ah 

+5{t" - t')Q{t"' - t')Q{t - t"')M2 

+S{t'" - t")Q{t" - t')Q{t - t")6{t" - t')M3 (4.43) 

+e(i"' - i")e(i - t"')e(i" - t')e{t - t")S{t" - t')M4 
+<d{t"' - t')Q{t" - t"')e{t" - t')e{t - t")S{t"' - t')Mz 

The coefficients Mi, • • • , Afs are presented in the appendix, (jA.28p . When this is inserted into p.lip . 
some of the time integrals can be executed. Moreover, we can apply an approximation since the time 



scale of the u-factors is much larger than that of the equilibrium time correlations. The upper limits of 
integrals over correlation functions can then be extended to infinity. The formula now reads: 

p3)o(i) - --p2/3ui(t)u2(i)u3(i) X 



/•OO nt nt — t 

x(/ dt'^loi23{t') + \im dt' dt"$2oi23(i',t")) (4-44) 

Jo *^°oJo Jo 



with $1 = Ml + M^ and $2 = M2 + M4 + M5. The detailed expression are given in the appendix, 
(|A.29|) . - This is the final form of the part of the third-order term which stems from the term S^^^ of the 
kernel function in the expression (j2.14p for the stress tensor. As has been mentioned earlyer, still to be 
added is the part of the second-order term which derives from p.3p . (j3.6p for e = 1. 



5 Relevant part of the second-order term 

When one treats the second-order term in the same way which led to p. lip , one obtains: 



00 



{D2h{t)^-\p^r^(^] ^ [dt'J^ dt"{K2)oi2(t,t',t")ui{t')uq{t") (5.1) 

In this section, the number index 2 contains a component index which has the value 1, if not stated 
otherwise. The pre-factor contains a part which results from switching to normalized variables in the 
correlation functions, uq is the Fourier transform of v? for the wave number variable k". The kernel 
function K2 is the sum of (|3.2p and p.Sp . taken at 6 = 60 for e = 1; (|3.10ap has been used to substitute 
the operator g: 

(i^2)oi2 (i, t', t") = -5{t" - t') {[g{t - t')f^]fl5hl) + 

+ e(t" - t')Q{t - t"){[g{t" - t')QCV6h1,g{t - t")ro]f*) (5.2) 

For the first term, we have: 

{[g{t - t')f,]fl5hl) = {[g{t - t')f^]flhl) - {[g{t ~ t')f,]fl){h*^) (5.3) 

The operator V in the second term in (|5.2p is evaluated: 



{[g{t"-t')QCV6h*^g{t~t")f^]rl) 

= -{[g{t t")f,]hlhl){[g{t" - t')h]fl) (5.4) 

The identity QCSh^ ~ —Qrs = —f^ has been used. Originally, the full form (|4.ip with 5h instead of h 
has been used in the first factor; but it is found that the terms with the constant parts (/12), {hz) vanish, 
so that the final form is (|5.4p . - For further preparation the formulas (|4.6p and (|4.7p . (|4.17p . (I4.2ip are 
used, integration over t" and an approximation because of different time scales is performed: 

{D2)o{t) = -\pip'^ {^^b^ ' ui{t)uq{t) j^ dt'{K2)oi2 (f) (5.5) 

K2 , now a function of t' only, is given in the appendix, (|A.30p . 



6 Result 

The quantity (15. 5p has to be added to (J4.44I) in order to obtain the part of the friction force which is 
third order in the velocity u. We denote it again D3: 
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P3)0(t) 



-p2/?(ui(t)u2(t)u3(i) X 



^CXD nt nt — t 

X( / dt'^loi23{t')+ \im dt' dt"$2oi23(t',0)) 







t— >oo 



+ (p/3apP|^)-5ui(t)7/(7(i) / dt'(i^2)012 (t')) 







(6.1) 



This formula is the final result of this paper. The next step is the determination of the equihbrium 
correlation functions C3 and C4. The MCT has been applied for this purpose. A formula for C3 is 
explicitly given in [5], while C4 has been calculated with the aid of the rules given in these papers. 
This will be described in a consecutive paper. - The transformation back to geometric space represented 
another major effort which resulted in closed formulas for the kernel functions. The formula obtained 
contains multiple space integrals; therefore, the hydrodynamic equation (I2.13P in this approximation is 
an integro-differential equation which is second-order in the local and third-order in the non-local part. 
This part of the investigation will be discussed in a another paper. 



A Appendix: Calculation of second order derivative 

In this appendix, formulas are presented which appear in the course of the analysis in the main part of 
the paper, but are too voluminous to be given there. - First, the terms of the second order functional 
derivative of the stress tensor kernel function S are written, as calculated from p.2p . (|3.5p . The detailed 
formula for the second derivative reads: 



S^Sabcdi^',t') 



(A 1) 

For shortness, on the left hand side, indices and variables are omitted. Second order derivatives with 
respect to the momentum part of the conjugated parameters are needed only (that is, the greek index e 
is switched to e). For each term, the first step is just the definition of the quantity. 



r^2^nfei) 



5h5h 



= -5{t" t')l3 tr{^^^Ml_[^(i',i)Q(i),^^(x)][Q(t')s,,(x')]5a,(x",i')} 



Sit'" - t')5[t" ~ t')P {[g{t', i)sac(x, t)] [Q{t')[Q{t')sM{^)]5a,{^\ i')]^a/(x"', t'))^,*' (A.2) 



5^S 
6b 6b 



-6{t" - t')p m', t)Q{i:)sU^)]Q{f)[Q{t')sM{^)] ^f^^'^l^'rj}^ )l.v 



= -6{t"' - t')6{t" - t')/3 {[g{t',t)Q{t)sac{^)]Q{t')sM{^'))LM{ae{^")6af{K"',t'))L,t' (A.3) 



^] ^ -6it" - t')l3 ([g(t^t)Q(t).,,(x)]Q(tO[ ^Jfff |.,) ^6.(xO]<5a,(x",tO)L,t> 

= -5{t'" -t')6{t" -t')x 

xP{[g{t',t)Q{t)saci^)]Q{t')[P{t')6afi^"',t')Qit')sM{^')]6a,{^",t'))L,t' (A.4) 



6^S 



(pA) 



= -6{t"-t')P{[g{t',t)Qit)Sac{^)] 



6Q{t') 



6bf{x"',t" 



[Q(i')sM(x')]<5ae(x",t'))L,*' 



6b 6b 

-6{t"' -t')6{t" -t')x 

X H[Q{t\m{t>ac{^)]r{t')6af{^'' ,t')Q{t')[Q{t')sM{^)]6a,{^\t'))L,t' (A.5) 
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5'^S 



n (p,5) 



5h5h 



= -5{t'' -t')(3 {[git' ^)^^j^^^sU^)]Q{t')[Q{t')sU^^ 

^ ~5{t"' ~ t)5{t" - t')x 

P {[g{t' ,t)V{t)5af{^" ,t)Q{t)sac{^)]Q{t')[Q{t')sM{^)]5a,{^' X))L,t' 
= 



(A.6) 



The last step is correct because of (|3.7p and QV — 



5^S 



>Sb 



''^^^""^'^'^^[^i^^^^*^'"'^''^^^^^'^^^^*'^'"'^'''^^^"'^^''"'*'^^^''' 



= -Sit" - t')e{t"' - t')e(t - t"')x 

x/3([e(i',t"')/:7'(t"')5a/(x"',t"')Q(i"')^(i"',t)Q(i).w(x)]Q(i')[Q(i')sM(x')]<5ae(x",t'))L.*' 



(A.7) 



rs^S.iM 



SbSb 



x/3tr{- 



^ Q{t" - t')e{t - t") X 



[g{t\t")cr{t")SaA,^",t")Q{t")g{t'\t)Q{t)sac{^)]Q{t')sM{^')} 



'%(x"',t"')' 

-(5(i"' - 1')©(^" - t')e{t - 1") X 

x/3([a(i',t")/:^(i")'^fle(x",i")Q(t")^(i",i)QWsac(x)][Q(t').'^bd(x'#a/(x"',i'))L,t' (A.8) 



r^2^n(/^2) 



56(56 



= e{t" - t')Q{t - t") X 



x/3([g(t^t'o/:p(t")fee(x^^")Q(t-ogft^t)Qft).,e(x)] ^^^gJ'|,„^ .,,(xO)L..> 

<5(i"' - t')e(i" - ^')©(^ - ^") X 

x/3(a[(<',t")/:^(i")'5«e(x",i")Q(t")a(i",i)Q(i)sac(x)]P(t')'5«/(x"',i')Q(^')s6d(x'))L,t'(A.9) 



r^2^n(/.3) 



(5656 



e(t" - t')Q{t - t") X 



<5Q(i) 



-Sac(x)]Q(i')s6d(x'))L,t' 



x/3([(?(i', t")£P(<")5ae(x", i") Q(i")^(^", ^^ (^,„^ ^,„) . 

5(i"' - t)<d{t" - t')Q{t - t") X 

x/3([^(i',t")/:P(t")^ae(x",i")Q(i")^(t",t)7'(i)5a/(x"',t)Q(t)sac(x)]Q(i')s6d(x'))L,t' 

(A.IO) 



5^S 



e{t" - t')Q{t - t") X 

<5Q(t") 



-^(t",t)2Wsac(x)]Q(i')sbd(x'))L,t' 



(56(56 
x/3([g(t',i")/:7'(i")<5ae(x",t")^^^(^,„^^„,). 

5(i"'-i")e(i"-i')©(^-i") X 

x/3(e[(<',i")/:7'(i")'5ae(x",t")7'(i")<5a/(x"',t")Q(i")e(i",OQWsac(x)]Q(t')sbd(x'))L,t'(A.ll) 
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r^2^n(/,5) 



Sb6b 

xi3{[g{t',t")cpit 



e{t" - t')<d{t - 1") X 



6/(x"',<"') 

-e(t" - t')(d{t - 1") X 
x/3([g(i',t")/:p(t"^'^^°'^'^"^^'^'*" 



Q(i")^(t",i)QWSac(x)]Q(t')s6d(x'))L,t' 



%(x"',t'' 

5(i"'-i")e(t"-t')e(t-t") X 

x/3([a(i',i")^7'(i")Q(i")^(i",t)QWsac(x)]Q(t')sM(x'))L.t'(ae(2:")'5a/(x"',i"))L,t" 

(A.12) 



5^S 



iffi) 



SbSb 

xp{[g{t\t")c 



= Q{t" - t')e{t - 1") X 

SV{t") 



Sa,{x",t")Q{t")g{t'\t)Q{t)sac{^)]Q{t')hd{^'))L,t' 



Sbf{x"',t"') 
-5{t"' - t")Q{t" - t')e{t - t") X 
xP{[g{t\t'')CV{t'')5af{^'\t'')Q{t'')5ae{x'\t'')Q{t'')g{t^^ 

X Q{t')sM{^))L,t' 



(A.13) 



5^S 



Q{t" - t')Q{t - t") X 



■Q(t)Sac(x)]Q(t')sM(x'))L,t' 



SbSb 

xp{[g{t',t")Cr{t")5a,{x",t")Q{t")^^ ^^,„^^,„^ 

Q{t"' - t")Q{t - t"')Q{t" - t')e{t - t") X 
xP{[g{t\t")CV{t")5a,{x" ,t")Q{t")x 

xg{t'' ,t''')CV{t''')5as{x''\t''')Q{t'''))g{t''' ,t)Q{t)sac{^)]Q{t')^^^^ (A.14) 



5'S 



Q{t" - t')Q{t - t") X 



xp{[ ^f^^^;f'^,^^^ ^^ 

Bit'" - t')e{t" - t"')e{t" - t')e{t - 1") x 

x/3([a(i', t'")CV{t"')Saf{K"', t"')Q{t"'))g{t"', t") x 

X CP{t'')5ae{x'' ,t'')Q{t'')g{t'' ,t)Q{t)saM)]Q{t')sM{-^'))L,r (A.15) 



Results (|A.10I) . (jA.12[) are found the same way as (IA.6p . - Below, the non-zero terms are presented 
in the form they assume after the reformulation in sec. 3: 



' 6'^R 



SbSb 



(p,i) 





SbSb 

S^R' 
SbSb 



^2j^,ipA) 



SbSb 



= -S{t"' - t')S{t" - t'){[Pro + Lg{t - t')h^][Qwlhl]hl)^ 

(P.2) 

= -S{t"' - t')S{t" - t') {[git - t')f^]fl)o{hlhl)^ 
= -S{t"' - t')S{t" - t') {[g{t - t')r^][Vflhl]hl)o 
= S{t'" - t')S{t" - t') {[Vr^ + Cg{t - t')ho]PhlQflh;)o 



(A.16) 
(A.17) 
(A.18) 
(A.19) 
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d'^R 



n (p,6) 



5b Sb 



= -Sit" - t')Q{t"' - t')Q{t - t"'){g[{t"' - t')QCVhlg{t - t"')f,]flhl)o 



(A.20) 



^2^. (/a) 



5b 5b 



= -Sit'" - t')Q{t" - t')e{t - t"){[Cg{t" - t')Vhlg{t - t")fo]f*/i*)o 



(A.21) 



^.^,im 



5b 5b 



5{t"' - t')Q{t" - t')Q{t - t"){[Cg{t" - t')rhlg{t - t")r^]Vflhl)^ 



(A.22) 



^.^.(/4) 



5b 5b 



5{t"' - t")Q{t" - t')Q{t - t"){[g{t" - t')QCPh;Phy{t - t")fo]fl)o 



(A.23) 



5b5b 



(/,6) 



= -5{t"' - t")Q(t" - t')Q{t - t"){[g(t" - t')QCVhlQhlg{t - t")ro]fl)^ 



(A.24) 



^.^.im 



5b 5b 



= e(t"' - t")Q{t - t"')Q{t" - t')Q{t - t") X 

X iW - t')QCVhlg{t"' - t")Q£rh*g{t - r)fo]f*}o (A.25) 



rs^ji.iM 



5b 5b 



^ e(i"' - t')e(t" - t"')e(i" - t')eit - 1") x 

X iW" - t')QCrhlg{t" - t"')QCVh;g{t - t")f„]fl)o (A.26) 



After evaluation of the projection operators in section 4, we obtain for the 11 terms, now cahed 
Ki, ■ ■ ■ , Kii (Symbol o omitted from ( )): 



(^i)o 

{K7)o 

iKs)o 






{[Cg{t - t,)ho]hlhl){rmhi) - {[Cg{t - ti)ho]rlh;hl) 

-{[g{t-t,)fon){h;hl) 

-{[git - ti)fo]hlhl){flhlh4) 

{[Cg{t - ti)ho]hl){flh*^hlhi) - {[Cg{t - ti)ho]hl){rlh*^h5){hlhihl) 

{[g{t-t3)foKhl){[g{ts-t,)h]flh;) 

-Mt - t2)fo]h;hl}{[Cgit2 - t,)hi]flhl) 

{[g{t - t2)f„]h*2K){[Cg{t2 - ti)h5]hl){flh;hi) 

-{[g{t - t2)fo]h;){[g{t2 - t,)h]n){h*2hl) 

+ {[git - t2)ro]hlhl}{[git2 - ti)r5]rl}{h*2hlhe} 

-{[Cg{t - t2)ho]hl){[g{t2 - ti)hn){hlhl) 

+ {[g{t - t2)f^]h*2hlhl){[g{t2 - t^)h]fl) 

+ {[g{t - t2)fo]h*2K){[g{t2 - ti)h]fl){hlKhe) 

{[g{t - t3)ro]hih;){[git2 - ti)hK){[git3 - t2)f7]h;hi) 

{[git - h)ro]h*2h*^){[g{t3 - ti)f5]fl){[g{t2 ~ t3)f7]hlht,) 



(A.27) 



In sec. 4, these terms are combined into 5 main parts: 
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(Ml) 0123 = -C237oi'5(i - t') 

-2Si24iii05iC3)534{t - t') + dt{C^)azi{t - t')) 
-Wi8((C'4)7823(i - t')Km + dtiCi)os23{t - t')) 

-l{K07dt{Ci)7123{t - t') + du(C4)0123{t - t')) 

(M2)0123 = ((C3)824(t " t'O^^OS + 9t(C3)o24(t - i")) 
x(k45(s513 - S135) 
+ 2c:;i6((C3)763(i" - t')k„ + at(C3)463(t" - <') 

+ «475t(C3)713(i" - i') + ^5,t(C3)413(i" " ^0) (A.28) 

(M3)oi23 = 75l(^08((C3)836(i-t")j625 + (C3)826(t-t")j635) 

+5t(C3)o26(i - i")j635 + 5t(C3)o36(i ^ t")j625 

-(C4)8235(i - i")«08 - at(C4)0235(t " *")) 
(M4)0123 - 75l(^08(C3)837(t-i"') + 5t(C3)037(i-i"')) 

x(^79(C3)925(i"' - t") + 9t(C3)725(t"' " t") 
(M5)0123 = 75l(«08(C3)827(t - i") + dt{C^)a27{t ~ t")) 

x(^79(C3)935(i" - n + at(C3)735(t" - <"') 

When the A/-quantities are inserted into (j3.1ip . some of the time integrations can be performed; 
finally, two parts of the kernel function remain: 

*l0123(i') = -5237oi'^(*') 

-2Si24(K05(C3)534(i') + 9* (C3)o34(t')) 
+2J62575l((«08((C3)836(i') + 5t(C3)o36(t')) 
-K5l((C4)8235(i')«08 + dt{C4)o23b{t')) 

-\{f^mdt{C^)7l23{t') + dtt{C4U23{t')) 

<^2oi23{t',t") = 275i(K08(C3)837(t')+at(C3)037(i')) (A.29) 

x(^79(C3)925(i") + 9t(C3)725(i")) 
+ ((C3)824(i')'^08+5,(C3)024(i')) 
x(k45(s513 - S135) 

+ 2cil6((C3)763(i")^47 + 5,(C3)463(t")) 
+ K47<9,(C3)713(i") + IdttiCsUisit")) 

In sec. 5, the part of the second-order (in b — 69) term which adds to the third-order velocity term is 
calculated. The kernel function reads: 



(i^2)012(i') - l0im{h2) 

+73l(5t(C3)032(0 + Ko4(C3)432(t)) 

+dtt{C3)012{t) + iKdt{C3)012\"2"=c{t)+i:Jl4dt{C3)0A2{t) 

+ K03{dt(C3)3i2{t)+iK(C3hi2\"2"=c{t)+<:Ju{C3)342{t)) (A.30) 
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